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CROSS  FIELD  JETTING  OF  ENERGETIC  IONS  PRODUCED 
BY  RAYLEIGH— TAYLOR  INSTABILITY 


I .  Introduction 

A  high  attitude  nuclear  explosion  (HANE)  can  lead  to  cross  field 
transport  of  debris  which  in  turn  produces  energetic  trapped  betas  on 
L  shells  considerably  higher  than  the  one  on  which  the  detonation  took 
place.  The  jetting  model  currently  used  in  the  SPECTER  code  is 
an  estimate  from  ad  hoc  arguments.  In  this  work  the  cross  field  transport 
of  charged  debris,  during  expansion  of  the  debris  sphere,  is  modeled 
by  considering  the  Rayleigh-Taylor  instability  in  a  linear  calculation. 

The  model  considers  three  phemonena.  Two  provide  energy  sources  for 
the  instability  and  the  third  provides  a  damping  mechanism.  The  damping 
mechanism  allows  estimates  of  minimum  wavelengths  for  the  modes  as  well 
as  critical  thickness  and  temperatures  for  onset  of  the  instability. 

The  first  phenomenom  which  can  initiate  the  Rayleigh-Taylor  instability 

is  the  laminar  like  accelerations  of  a  percentage  of  the  ions  within  the 

shock1 .  Essentially  this  leads  to  a  three  species  description  of  the 

plasma.  Particle  simulations  of  an  advancing  shock  have  demonstrated  the 

existence  of  these  accelerated  ions  which  have  densities  of  approximately 

10#  of  the  ion  density  in  the  shock.  The  third  species  causes  a  space 

charge  with  the  resultant  electric  field.  The  result  of  this  is  the 

Rayleigh-Taylor  instability.  The  second  source  of  energy  for  the  instability 

is  the  curvature  of  the  magnetic  field  lines  caused  by  the  bubble  expansion. 

This  force  is  more  apparent  during  early  phases  of  the  bubble  evolution. 

Finally,  there  is  damping  of  the  instability  provided  by  finite  Larmor  radius 

corrections  to  the  momentum  equations  (i.e.  magnetic  viscosity).  The  damping  is  a 
Note:  Manuscript  submitted  June  25,  1979. 
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function  of  the  instability  wavelength  and  is  directly  coupled  to  the 

thickness  of  the  debris  shell  as  well  as  the  perpendicular  pressure  of  the  ctebris. 

In  this  report  a  simplified  model  for  prediction  of  these  effects 
will  be  presented.  Details  of  the  calculation  of  all  three  phenomena  are 
discussed  in  Appendix  A  and  B  and  some  sample  estimate  of  growth  rates, 
wavelengths  and  mass  transport  are  presented  in  Appendix  C. 

II .  Model 

In  this  model  we  make  use  of  the  growth  and  damping  rates  represented 

by  Eq .  A-21  and  Eqs.  B-9,  D,  11.  Use  is  also  made  of  some  scaling  relations  provided 
by  Dr.  Robert  Clark.  These  relations  have  been  previously  provided  under 
this  contract.  As  in  Appendix  B,  we  work  in  cylindrical  coordinates  where 

*  A 

the  B  field  is  in  the  9  direction  and  displacements  are  in  the  r  and  z 


direction,  Fig.  1. 

The  growth  rate  from  the  laminar  acceleration  is  represented  by  the 
following  equations 
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' L  z  A 
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where  , |—  |  ,  Lm  =  |g— |  ,  is  the  Alfven  spaed  and  f  is  the 

o  Po  o 

fraction  of  ions  accelerated  in  the  shock,  typically  ICfjh.  In  this  report 


L  is  taken  equal  to  L  .  Equation  II- 1  provides  one  of  the  restrictions  on 
n  m 


Xz» 


X  <2ttL  • 
z  n 


II-2 


It  should  also  be  noted  that  as  X_  -*  o,y  -»  constant.  The  damping  rate 

Z  L 
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produced  by  finite  larmor  radius  effects  is  calculated  in  Appendix B. 
and  is  found  to  be 


o  pK  k  k 
x  r  z 


P  U)  . 
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H-5 


The  growth  rate  due  to  the  centrifugal  acceleration,  calculated  in 
Appendix  B,  is  represented  by  one  of  two  relations: 
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Equations  11-1,3,5  lead  to  estimate  of  the  minimum  wavelength  in  the  z 
direction  ie:  the  condition  for  growth  is 


yL  +  yc  >  V  11-6 

Where  Eq.  II-5  is  the  limit  most  likely  to  occur  in  this  problem, 
using  equations  II-i,  II-3>  and  II-5>  one  obtains  the  following  relation 


or 
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Equation  II-2  and  II-8  provide  the  upper  and  lower  bound  on  \z>  One 
additional  condition  is  required  to  complete  the  discription  that  is 
\z  <  tt  R  .  This  implies  that  the  wavelength  must  be  less  than  the  upper 
half  of  the  circumference  of  the  bubble.  These  three  conditions  determine 
the  time  of  onset  for  the  instability  and  provide  additional  information 
as  to  possible  orientation  of  region  of  maximum  instability. 

The  next  estimate  required  is  the  amount  of  mass  involved  in  the  instability. 
If  we  consider  a  typical  sinusoidal  oscillation  as  the  perturbation,  an 
estimate  of  amount  of  debris  can  be  achieved.  See  figure  2. 


If  £r  C  X  A  and  \/k  <  R  we  can  assume  angles  0  and  ^  small  and  arrive 
at  the  following  relation 


<r 


II-9 


where  ntust  be  less  than  A,  the  coupling  shell  thickness.  In  order 

for  the  instability  to  remain  in  the  linear  regime  it  has  been  shown2  that 

C  <  -2X,  .  We  will  take  C  ~  .1  X  .  As  the  instability  progresses  into 
r  z  r  z 

the  nonlinear  regime  most  of  the  mass  involved  can  be  reasonably  es¬ 
timated  by  the  area  within  the  linear  regime.  The  area  subtended  by 
a  triangle  which  covers  the  ejected  mass  regions  of  a  single  oscillation 
is  approximately 

A~^zcr~To  11-10 


The  total  mass  ejected  per  lobe  is  then  estimated  by  poAigand  the  total 
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number  of  debris  particles  is  estimated  by 


°T  ~  W 


II-ll 


The  scale  length,  ,  is  considered  to  be  the  wavelength  in  the  0 
direction.  Since  the  Rayleigh-Taylor  instability  is  a  flute  type  mode 

typically  k*  B  w  0,  however  there  can  be  finite  k.,  or  in  this  case  k 

i  6 

as  long  as  kjj/k^  <  (me/m^)  .  This  establishes  an  upper  bound  on  .  A 

A 

lower  bound  is  approximately  half  the  hemispheric  distance  in  the  0  direction 

so  the  minimum  =tt*  Here  we  will  take  the  conservative  estimate 
£q  =  .  The  volume  of  debris  involved  in  each  outward  jetting 

lobe  can  now  be  calculated  from  Eq.  II-ll. 

The  scaling  used  in  this  model  was  obtained  from  numerical  calcu¬ 
lations  by  Dr.  Robert  Clark.  We  begin  with  the  usual  kinetic  energy 
re  la t ion 


Yk ' 1/2  W 


11-12 


where  is  the  kinetic  yield  in  ergs. 
Assuming  complete  coupling  and  p^(R)  ~  const 


Mu  =  [M  ttR3pa^u 


more  exactly  Mu  =  (Mq  +  \tt  fR  R'e  pA(R’)dR')u 
*°  °  °  **0 

If  we  define  R  3  =  M  /(4/3  rrp.  )  as  the  radius  at  which  a  weapon  mass  of 

O  A 

air  is  swept  out  then  the  following  scaling  appears: 


u/u  -  1/(1  +  R3/R*3) 


11-13 
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|~  =  l  +  r3/r*3 

o 


11-14 


UQt  =  R(1  +  1/4  R  %  3) 


11-15 


Making  use  of  the  marginal  stability  criterion  for  the  ion-ion  instability 
and  assuming  is  less  than  unity,  a  scaling  relatim  can  be  derived  for  the 

B  field  in  the  shock  as  a  function  of  radius 

u(^TTp. 


B 


o  2.5 

Using  Eq.  11-15  we  obtain  the  relation 


B  ~  <*  u  (Imp  )*  /2. 5(  1  +  R3/R*3 ) 
o  o  A 


11-16 


where  BQ  2:  B^  must  hold.  This  relation  is  an  underestimate  of  the 
B  field  as  the  ion-ion  instability  has  already  begun.  Comparisons  with 
simulations  show  that  a  factor  of  3  is  necessary  to  bring  this  relation 
into  agreement  with  the  simulations.  Therefore,  the  factor  a  is  set 
equal  to  3  in  recognition  of  this  fact. 

The  scaling  of  temperature  is  obtained  by  noting  that  T^  <  T^inetic 

must  hold  because  Landau  damping  for  the  ion-ion  instability  occurs 


where  V 


fW. 

V.  <  V.  s  W — - 
A  i  ^  m. 


We  then  use  the  following  scaling: 


mu  ^kinetic 

(kTi/m1u2)  ~  .2 


therefore  kT^  ~  .04  m^u2  ~  .04  kTkinetic  11-17 

Substituting  Eqs.  II-16,  11-17,  11-13  into  Eq.  II-8  we  obtain 


6 


12n2k  L  ( .04)m, u2p. 
r  n'  i  yA 


u)  .  B 
ci  o 


(kTT0.  )$ 


or 


5nLnkr(.04)(2.5)2 

Z  «2 ( ^PA  ^°a ^ po  ^  ^  +  d 

This  minimum  wavelength  can  now  be  estimated  if  k^  and  are  known. 

We  estimate  the  maximum  X to  be  the  coupling  shell  thickness,  A, and 

estimate  the  density  and  magnetic  scale  lengths  to  be  the  e-folding 

distance  of  half  the  coupling  shell  thickness,  Lr  ~  .37  A/2.  From 

simulations  c£  Dr.  Clark  it  is  found  that  the  coupling  shell  thickness 

with  respect  to  the  magnetic  bubble  radius  varies  with  time.  For  reasons 

to  be  discussed  in  the  following  paragraph  the  instability  will  have  its 

onset  when  approximately  one  weapon  mass  has  been  swept  out.  For  Starfish 

* 

the  time  at  which  R  =  R  is  approximately  120  msec,  at  this  time  A  ~  .025R. 

* 

The  estimate  that  the  instability  onset  begins  at  approximately  R=R 
is  determined  by  Eq .  II-2.  At  this  point  in  the  dynamics  a  variety  of 
effects  appear  but  the  dominant  effects  are  the  drop  in  temperature  and 
continued  increase  in  the  coupling  shell  thickness  allowing  the  condition  on 
Eq.  II-l  to  be  met.  Applying  this  scaling  for  kr  and  Ln,Eq.  II-2  and 
Eq.  II- 18  can  be  written  as 


m.  c 

f 


11-18 


Xz  <  0.37  A/2 


11-19 
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and 


x  >  6tt2(.37)(.04)(2.$)2 

Z 

An  additional  piece  of  information  is  also  available  from  Eq.  II- 19. 

The  coupling  shell  thickness  is  a  function  of  angle  with  respect  to  the 
magnetic  field  approximated  by 

A(R,9  )  =  A(R)/sin9,  0  >  0O 


where  9q  is  the  angle  at  which  coupling  fails,  9q  ~/A/R. 

Equation  II- 19  provides  the  major  restriction  to  the  onset  time  of  the  in¬ 
stability.  As  can  be  seen,  increasing  A  leads  to  earlier  onset  of  the  mode, 
therefore  there  will  be  a  preferred  direction  to  the  largest  modes.  Since 
the  field  lines  are  frozen  into  the  plasma  and  the  instability  is  a  flute 
instability  k*  B  =  0,  the  region  of  max.  A,  where  the  expansion  is  per¬ 
pendicular  to  the  frozen  field  line,  v'ill  be  some  angle  between  field 
alignment,  9=0  and  perpendicular  to  the  ambient  field  line  9  =  tt/2,  but 
will  easily  satisfy  the  condition  9  >  9.  For  symmetrical 

expansion  into  a  B  field  in  a  constant  density  atmosphere,  one  expects  two 
regions  of  earliest  onset,  see  Fig.  3.  However,  the  atmosphere  is  exponential 
and  therefore  one  has  a  preferred  direction  along  the  ambient  density  gradient 
where  A  will  actually  be  the  largest.  Fig.  4.  This  argument  suggests  that  the 
region  of  maximum  growth  will  be  in  a  direction  slightly  southward  and 
essentially  upward  for  Starfish,  Fig.  4.  This  model  predicts  onset  of  the 
Rayleigh-Taylor  instability  at  a  time  when  R  =  R  .  At  this  time  the  wavelength 
is  estimated  to  be 


x  _ hhQ _ 

z  e  o2(4ttpa)£  [(pa/pq)*  +  (f/Z)i 


11-21 
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and  the  number  of  debris  atoms  injected  is  represented  by 


ND  " 


20 


11-22 


Here  n^  is  taken  to  be  0.5  (ng)  with  ng  ~  4  r.^  where  n^  is  the  ambient 
ion  density,  R  is  as  usual  the  radius  o'f  the  bubble  and  0 ^  is  the  ambient 
mass  density  of  the  ions.  Generally  all  quantities  labeled  /  are  defined 
in  the  coupling  shell  and  the  ambient  is  defined  to  be  the  air  located  at 
the  shock. 

III.  Summary 

The  Rayleigh-Taylor  instability  is  found  to  be  unstable  with  growth 
rates  determined  by  Eqs.  11-1,3,5  and  wavelengths  determined  by  II-2 
and  11-21.  The  wavelength  and  the  total  number  of  debris  particles  trans¬ 
ported  across  the  field  lines,  NQ,  are  found  to  be  functions  of  altitude. 

The  instability  appears  to  onset  When  a  weapon  mass  of  air  has  been  swept 
out  by  the  expanding  shock  front.  From  considerations  of  the  relative 
thickness,  A,  of  the  coupling  shell  a  directional  dependence  can  be  surmised. 
The  direction  most  optimal  for  the  instability  appears  to  be  oriented  in 
the  direction  of  the  ambient  density  gradient. 

For  the  purpose  of  modeling  the  following  relations  appear  to  produce  the 
appropriate  results.  The  wavelength  is 


5-^8  m^c 

Xz  ^(^Pa^Pa/Po^  +  (//2)i] 

where  a  =  3,  p  is  the  mass  density  in  the  coupling  shell,  /  =  0.1,p.  is  the 

O  A 

ambient  density  at  the  coupling  shell  and  m  =  27  m  .  The  total  number  of  debris 

1  0 
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particles  transported  across  the  field  lines  is  computed  by  the  relation 


ND  = 


20 


V 


with  n_  =  0.5  n  where  n  is  the  density  in  the  coupling  shell  and  n  ~  4  n 
D  S  S  S  A 

where  n  is  the  ambient  mass  density  at  the  coupling  shell  altitude.  The 
A 

„  #  * 

9  scale  length  is  taken  to  be  ig  =  tt  R/2  where  R  R  and  R  is  the  radius 

at  which  one  weapon  mass  of  air  has  been  swept  out. 
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Appendix  A.  Linear  Calculation  of  Rayleigh-Taylor  Instability 


Driven  by  Shock  Acceleration 

In  this  calculation  it  is  assumed  that  the  shock  is  moving 
perpendicular  to  the  magnetic  field  lines.  The  shock,  spatial  gradients 
and  force  are  taken  in  the  2  direction  with  the  magnetic  field  in  the  £ 
direction.  Three  species  are  considered.  The  electrons  and  ions 
moving  in  the  shock  frame  of  reference  and  a  third  species  com¬ 
prising  approx.  10$  of  the  density,  accelerated  within  the  shock  with 

,  m2V 

a  force  F  =  £  — -  ,  where  V  is  the  Alfven  speed  and,  L  is  the 

L  «.  m 

m 

magnetic  field  scale  length.  Further,  it  is  assumed  that  the  magnetic 
field  is  frozen  into  the  plasma. 

Five  variables  are  unknown,  the  perturbed  velocities  for  the  two 
ion  species  and  the  three  perturbed  densities.  Therefore,  five 
equations  are  necessary.  Consider  the  general  momentum  equation: 


dv .  n . 

n.m.  +  n.m.(v.  -V)v.  =  q  .  — ^  v ,  x  B  +  n.q.E  -  VP  +  n.F.  A-l 

J  J  St  j  j'-j  -'-j  c  -j  -  pi-  -j  j  j 


and  apply  the  following  scaling  cu.^»  ,V  •  V  and  u)ce  »  This 

implies  that  inertial  terms  can  be  neglected  and  thatve^/u)ce  «  v^e/u>C£ 


represents  collisional  drag  or  acceleration  depending  on  the  species 


and  is  in  the  z  direction.  With  this  ordering  the  electron  momentum 


Eq.  becomes: 
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The  ion  momentum  equation  to  lowest  order  in  the  cyclotron  frequency 
becomes: 


n. 


qi  ~(-i  x  ?)+  Vi-  -  +  niFi  with  i  =  1,2 


A-3 


so  that  the  lowest  order  in  velocity  is: 


/  \  B  x  VP.  Bx  F. 

V  =  +  c  1 - li  _  c  _ li 

qiniB2  qiB2 


A-4 


To  the  next  order  one  obtains 

(o) 


nimi  “Tt“  +  (^i°-  ^Yi(o)]  =  Vi  [?  +  ~T~  ] "  ^i+V-i 


A-5 


In  addition  to  Eqs.  (A2  and  5)  we  use  the  continuity  equation 

3n. 

^+V.(njvj)=0  a-6 

and  the  ideal  gas  assumption 

P .  =  n . T .  A-7 

J  J  J  A  ' 

We  shall  assume  charge  neutrality,  V-  J  =0,  and  incompressibility, 

Y  *  Yj  =  °*  Further,  we  substitute  v.(o^  for  v.  in  the  equations  and 
then  write  v^  to  represent  v^  \  The  assumption  that  °^is  equal 
to  the  is  valid  if  cB  x  E  cB  x  VPe  cB  x  F  which  is 

- —  »  — : — 11  "  _i - e 

n  q  B2  q  B2 

e^e  Me 

generally  true.  Therefore,  Equation  (A-6)  becomes 


Vnj 


0 


A-8 
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Adding  Eqs.  A-2  and  5  we  obtain 
dv  i  dv 


av  i  uvp  \ 

ni"i  jr  +  V*  ar  -  (».s  +  "i’i  +  V*’5  +  (n.V«  +  “iVii'w*)*  ~ 

-  ’(p.  +  Pj  +  p2)  +  n/.„  +  ni!i  +  aJs 


Recall  that  7*J  =  0  implies  that  one  can  write 

J  =  -e  x  Vf  A -10 

where  i>  is  a  scaler  potential  function  such  that  J  x  B  =  -BV^. 

Using  Eq.  (A-10 )  in  Eq .  (A-9 )  and  making  useof  quasineutrality,  we  obtain 

dv!  dv2 

n.m  TT~  +  nomo  TZ~  =  +  £.in ■  f4  A- 11 

11  dt  22  dt  j  J  j 

A 

where  P  =  Pe  +  P^  +  £2  +  5^  and  we  have  assumed  this  problem  is  electro¬ 
static.  At  this  point  two  ways  to  approach  the  problem  are  available 
either  linearize  Eq.  A-TI  as  it  stands  using  =  f  and 

A  A 

keep  fourier  components  in  both  directions,  x  and  z,  or  move  to  a  center 
of  mass  of  motion  for  the  ion  species  and  then  linearize.  If  the  first 
option  is  applied  one  must  use  energy  conservation  to  close  the  set  of 
equations.  The  result,  to  first  order  in  k,  is  a  complex  cubic  equation 
with  complicated  coefficients.  The  second  method  which  will  be  followed 
here  requires  the  following  definitions: 


e"jVj 


where  p  =  Zn^m^ 


F  =  ^jnjFj  and  the  perturbation  takes  the 
form  f(x  z)  =  F(z)ei^xX  ut)  # 
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I  w 


Incompressibility  is  still  required  for  the  center  of  mass  velocity  as 
well  as  the  individual  ion  velocities.  Equation  (A-ll)  becomes 


P 


dV 

__  _  Vp 


+  F 


(A-12) 


when  V*V  =  7*V^  -  =  0  is  used.  Taking  the  curl  of  Eq.  (A-12),  and 

noting  that  Vq  is  in  only  the  x  direction  since  the  calculation  is  done 

A 

in  the  reference  frame  of  shock  moving  in  z,  and  linearizing  the  resulting 
equation  leads  to 


dv  -I  -  dV 

-  ik  1  +  - - 

x  dt  3z  dt 


dV. 


lx 


dt 


-  ik  F 
x-z 


(A-13) 


Applying  V*V  =  0  to  equation  (A-13)  leads  to  a  second  order  equation 


or 


,  3V. 

— —  p  Z  -  k2p  Vn 

3z  o  3z  x  o  lz 


ik  F. 
x-lz 

(oj-k  V  ) 
x  ox 


(A-14) 


Since  we  are  in  the  shock  frame  all  accelerations  on  the  shock  are  zero 

with  the  exception  of  those  generated  within  the  shock  itself.  The  force 

on  the  second  ion  species  is  due  to  a  type  of  laminar  acceleration  within 

the  shock. ^  Therefore,  the  force  term  reduces  to  F  ,  >  n„,F„.  In  order 

zl  z  12 

to  proceed  the  force  term  is  multiplied  and  divided  by  a  factor  of  p,  where 
the  ratio  of  the  perturbed  density  to  the  perturbed  center  of  mass  density 
is  taken  as  zeroth  order  so  that 


zl 


n21F2 


n21F2 


P  F 


2* 


Linearizing  the  center  of  mass  continuity  equation 

+  V*Vp  ”  0 


(A-15) 
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one  obtains  -i(w-k  Vox)p,  +  V,  p'  *  0.  Solving  for  p.  and  substituting 
x  1  lz  o  1 

into  the  force  term  one  can  rewrite  Eq.  (A-14) 


3  p«4t  V,  -  k2  v,_  + 


k2F  P'V. 

X  o  o  lz 


3z  °  3z  lz  Kx  o  lz  (oH-k  V  )2 

x  ox 


0. 


(A-16) 


The  solution  to  this  equation  is  easily  obtain  and  has  the  form. 


-  o/2p  ±iKz 

V,  =  e  0  e 

lz 


where 


„  F  k  p*  , 

/  ,  \2  ox  O/p 

M  +  ' k*  • 


K* 


The  boundary  conditions  applied  to  the  oscillatory  portion  of  the  solution, 
±iKz 

q  =  e  ,  leads  to  the  quantization  condition. 


K2d2  =  n2TT2 


(A-17) 


where  the  boundary  conditions  are  q  =  0  at  z  =  0  and  q  =  0  and  z  =  d.  With 
z  measured  from  the  peak  of  the  shock  profile  and  d  is  the  distance  to  the 
foot  of  the  shock  from  z  ”  0.  Using  Eq.  (A-17)  leads  directly  to  the 
equation 


'  F  p 


(“  -  Vox)  -  kx 


o  o/p 


1/2 


(A-18) 


k2 

X  , 


Now  p"/p  “  -1/i  on  the  front  side  of  the  shock  and  there  will  be  a  growth 
o  o 

rate  from  this  mechanism  if 


k2  + 


$ »  M  -  ^ 


(A-19) 
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If  (^p)2  >  ^2  the  mode  wi:L1  become  convective.  One  can  then  estimate  the 
order  of  n  that  will  allow  the  mode  to  grow.  Since  SL  =  de  \  then  mr  <  -|y 
must  be  met.  This  means  only  the  n  =  0  terms  will  allow  growing  modes. 
Therefore,  the  wavelength  of  the  instability  must  obey 

X  <  2n«.  .  (A-2C) 

x 

Providing  this  condition  is  met  the  growth  rate  becomes 

Pq  |  ^ 

2L"(k5-(cy'o)r) 

where  f  is  the  fraction  of  accelerated  ions  and  the  scale  length  of 
the  magnetic  field.  In  the  cylindrical  coordinates  employed  in  the  main 
body  of  this  work  kx-*kz. 


,1/2 


(A-21) 
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Appendix  B.  Rayleigh-Taylor  Instability  Driven  by  Centrifugal  Forces 

To  calculate  the  effects  of  centrifugal  forces  and  the  Larmor  radius  effect, 
it  is  necessary  to  work  in  cylindrical  or  spherical  coordinates.  For  this 
problem  cylindrical  coordinates  allows  the  necessary  physics  to  be  described; 
therefore,  we  will  work  in  this  system.  Here  we  will  work  in  Lagrangian 
with  the  various  quantities  described  as  follows: 


(0 ,Bg (R  ),0  )  R(t)  =  R(to)  =  C(t) 


B 

-o 

R  (t  ) 

u  (t  )  =  -f-2 

-o  o  dt 


d£(t) 

=  0  ^(0  =  —  =  1«C 


where  £  =  (gr,  0,  Cz)  and  f1  =  ■tot^ 

Four  equations  are  necessary,  we  use:  Faradays  law 
1 

-  —  =  V  x  u  x  B  ,  where  Ohms  Law 

C  ot  —  —  — 

E+uxB/c  =0  has  been  used,  the  momentum  equation 


du 

0  at  = 


(=  +  |?)  +  w  }  + 


u  x  B 


B-l 


B-2 


,  the  continuity  equation 

dn  +  v  .(n  u)  =  0  where  we  assume  V  .u  =  o 
dt 

and  finally  we  assume  for  the  pressure  equation,  the  ideal  gas  law 

P  =  nT. 


B-5 


B-4 


Equations  B-l, 4  are  now  linearized.  The  off-diagonal  pressure  elements 
provide  the  finite  Larmor  radius  effects.  The  off-diagonal  pressure  are 
represented  in  the  following  fashion, 

-du  3u 


_  r VaV| 

TT  -  IT  ~  -  T - 1-  - - I  , 

rr  zz  a>ct  I  dr  dz  J 


A?  -  £]  . 


B-5 


TT  ~  P 
rz 
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where  tt ^  =  TTrz .  The  set  of  linearized  equations  consisting  of  two  com¬ 
ponents  of  the  momentum  equation,  Faradays  law  and  the  continuity  equation 
are  the  placed  in  matrix  form  and  the  determinate  set  to  zero  to  determine 
the  frequency.  After  considerable  algebra  a  dispersion  relation  accurate 
to  the  second  order  in  k  is  obtained 


G  ±  (g2-4a)£ 


k  P  ,*  /  ,  an  \ 

=  -5^-  (3/r  +  -i-  — 2  -  (3/2  k  )) 
®ci  V  no  3r  r  / 

P  *k  /  -  an  \ 

h  U/T  +  i_  °) 

u,ci  V  no  3r  / 


«  s  1  Txkr 


‘ 3  •  kl(1/T  *&)**-  ik's) 

K  =  P  72(1)  .  /k  fl/r  +  1/r  -k  k  \ 

i  ci  l  z  y  o.3r  f  r  zl 


Neglecting  terms  of  order  1/r2  and  k3 


?  *  k  / 

3  U)  .  I-' 

o  ci  \ 


G  ~  — — — —  3/r  +  1/r  ^  -  i  3/2  k 

0  w  .  V  o  3r  "  i 
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Equation  B-6  is  in  a  form  to  be  handled  numerically  since  all  of  the  components 
are  either  zeroth  order  terms  or  must  obey  scaling  relations  which  have  to  be 
postulated  to  proceed  any  further. 

The  functional  dependence  of  the  growth  rate  for  the  wave  can  be  ob¬ 
tained  by  approximating  or  expanding  some  of  the  components  of  Eq.  B-l6. 

For  example,  if  one  assumes  the  gyro-radius  terms  are  smaller  than  the  curva¬ 


ture  term  then 


>G=or(G2-4A)i-12Ai(1-GS/8A) 


-  !  ±  iA*  (1- 


(1-G2/8A)  . 


The  variable.  A,  can  be  approximated  in  three  separate  forms: 


case  1 


case  2 


[>-=h:  ((&)'*&)')] 

*-  r  x  o 


1  i  B 

^~±  - ^ 

2r(pon)‘ 


(  Bo/r) 


case  3 


k  ~  *  3r 


Since  G  is  complex  there  will  be  both  a  real  frequency  and  a  growth  rate. 
The  three  cases  lead  to  the  following  solutions: 


case  1  linA  >  R-eA 


y  _  _  i  I^_  k  k  t  /Vi _l_  r  2  +  (\)]\ 

Po^i  r  2  V2mi  no  M  8nk  T  n'  L  0  Vr/J/ 

r  A  o  B-c 


case  2 


As  can  be  seen  from  these  relations,  the  finite  larmor  radius  terms 
provide  a  damping  term  and  the  curvature  affects  provide  a  growing  term. 

For  the  work  of  current  interest  only  case  1  or  5  are  of  interest.  Typically 
case  3  appears  to  be  the  dominate  limit. 
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Appendix  C 
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8  7 

then  l  1.09  x  10  .  If  we  use  &  *>  ttR  then  1  'v  6.97  x  10.  Actually 

0  0  0 

because  of  the  effects  of  the  exponential  atmosphere  in  A  as  discussed 

in  section  II.  We  estimate  £  'v  ttR/2  and  use  this  as  a  lower  bound  on 

u 

Then  ^  3.5  x  102  and 

U  u 

Nd  'v  5.3  x  1025  a  =  3  C-3 

min 

ND  'v  1.7  x  102^  a  =  3  C-4 

max 


using  the  upper  bound  as  £  . 

0 

Using  the  result  from  C-l  the  total  number  of  wavelengths  that  will 
fit  on  the  bubble  is 


Nw  =  p  'v  2.23  x  102  ,  a  =  3 

05 

N  *  25,  a  =  1 
w 


The  growth  and  damping  rate  have  the  following  magnitudes.  First  we 
calculate  some  useful  parameters 

VA  'v  .  3/ (4ttpa)  1  /2  %  3.2  x  107 

R  'v  R* 

u  'v  u  /2 
o 

2 Y  \  1/2 

-r)  '  !-7  x  108 

B  'v  au„  (4irp.  )1/2/5 
o  o  a 
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From  Eq.  II-l  and  5 


yi  +  ycn 


it  [( A/p.)  +(£/2)  /J-  ui  * 


From  Eq.  I 1-3 

Yn  *  -  1.5  (.04)  u2  ^ 

D  A  Az  Be 

YD  ^  -  1.21  x  102 

Y  =  YD  +  YL  +  Yc  ^  7.1  x  101 


06 


This  mode  will  e-fold  quite  rapidly  after  onset  of  the  instability. 
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Fig.  2  -  Linear  deformation  of  the  plasma  surface  during  growth  of 
Rayleigh-Taylor  instability.  |r  is  the  linear  displacement,  typically 
less  than  20%  of  X. 
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01CY  ATTN  DOC  CON  FOR  A.  LA80WITZ 

EG&G,  INC. 

LOS  ALAMOS  division 

p.  O.  BOX  809 

los  alamos,  nm  85544 

01CY  ATTN  DOC  CON  FOR  j.  BREEDLOVE 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

los  alamos  scientific  laboratory 


p.  0.  BOX  1663 
LOS  ALAMOS,  nm  87545 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 

SAND  I A  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  NM  87115 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 
01CY  ATTN  DOC  CON 


FOR  J.  WOLCOTT 
FOR  R.  F.  TASCHEK 
FOR  E.  JONES 
FOR  J.  .MALIK 
FOR  R.  JEFFRIES 
FOR  J.  ZINN 
FOR  P.  KEATON 
FOR  D.  wESTERVELT 


FOR  J.  .MARTIN 
FOR  W.  BROWN 
FOR  A.  ThCRNBROUGh 
FOR  T.  WRIGHT 
FOR  D.  OAHLGREN 
FOR  3141 

FOR  SPACE  PROuECT  DIV 
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SANOIA  LABORATORIES 
LIVERMORE  MORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

OVCT  ATTN  DOC  CON  FOR  d.  moRP-Cy 
OICY  ATTN  OOC  CON  POR  T.  COCK 

Of PICE  OF  military  application 

department  of  ENERGY 

-ASnlN&TQN,  O.C.  20545 

01CY  ATTN  >X  CON  FOR  3.  GAlE 


OTmER  GOVERNMENT 


CENTRA*.  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  Rm  5G48,  HQ  BlDG 
-ASmIN&tqn,  O.C.  20505 

OICT  ATTN  OS1/PSID  RM  5F  19 

department  of  commerce 
NATJONAc  3uREAl  OF  STANDARDS 
WASHINGTON,  O.C .  20254 

CALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  MOORE 

DEPART, MENT  OF  TRANSPORTATION 
OFFICE  OF  The  SECRETARY 
tad-44.1,  room  10402-0 
400  7Th  STREET,  S.w. 

WASHINGTON,  D.C.  20590 
01CY  ATTN  R.  LEWIS 
01CY  ATTN  R.  DOHERTY 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  S  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  i  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 

01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  Or  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 

P.  0.  BOX  92957 

LOS  ANGclES,  CA  90009 

01CY  ATTN  I  .  GARFuf*EL 
01CY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOER 
01CY  ATTN  N.  STOCKWELL 
01CY  ATTN  D.  OLSEN 
01CY  ATTN  J.  CARTER 
01CY  ATTN  F.  MORSE 
OiCY  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01603 


OICY  ATTN  RADIO  SCIENCES 

OICY 

ATTN 

OICY 

ATTN 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 

OICY 

ATTN 

P.  0.  BOX  983 

OICY 

ATTN 

BERKELEY,  CA  94701 

OICY 

ATTN 

OICY  ATTN  J.  WORKMAN 

OiCY 

ATTN 

OICY 

ATTN 

OICY 

ATTN 

BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KE^EY 

CALIFORNIA  AT  SAN  DIEGO,  UNlV  OF 
IPAPS,  B-019 
LA  JOLLA,  CA  92093 

OICY  ATTN  HENRY  G.  BOOKER 

3RON  ENGINEERING  COMPANY,  INC. 

CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  A*.  35807 

OICY  ATTN  ROMEO  A.  DELIBER1S 

Charles  star.<  draper  laboratory,  inc. 

555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  ATT-N  D.  6.  COX 
OICY  ATTN  U.  P.  GILMORE 

COMPUTES  SCIENCES  CORPORATION 
6565  ARLINGTON  3LVD 
falls  Church,  va  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 
LI NTh I  CUV  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

department  of  electrical  engineering 

ITHACA,  NY  14850 

OICY  ATTN  D.  T.  FARLEY  JR 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75060 

OICY  ATTN  H.  LOGSTON 

OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94066 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
QiCY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  &  COMMUNICATIONS  CORP 
3939  FABIAN  wAY 

palo  alto,  ca  94303 

OICY  ATTN  d.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICT  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  F.  RE  I BERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
316  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
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COMMANDER 

MARRY  OIAMONO  LABORATORIES 
department  of  the  army 
2800  POWDER  MILL  ROAD 
AOELPHI,  MD  20783 

(CNWDI-INNER  ENVELOPE:  ATTN:  DEL^D-RBh) 
01CY  ATTN  DELMD-TI  M.  WEINER 

01CY  ATTN  DELMD-RB  R.  WILLIAMS 

01CY  ATTN  DELHD-NP  F.  WIMENJTZ 

01CY  ATTN  DELMO-NP  C.  MQAZED 


COMMAMOER 

U.S.  ARMY  COMM-ELEC  ENGRG  INSTAL  agy 
FT.  HUACHUCA,  AZ  85613 

01CY  ATTN  CCC-EMEO  GEORGE  LANE 


COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  S  TECH  CTR 
220  7Th  STREET,  NE 
CHARLOTTESVILLE,  VA  22901 
01CY  ATTN  DRXST-SD 
01CY  ATTN  R.  JONES 


CGMMANDER 

naval  surface  weapons  center 

JAmcvjREN  laboratory 
OAhlGREN,  vA  22448 

01CY  ATTN  CODE  DF-14  R.  SUTLER 

COMMANDING  UFFi C5« 

NAvY  SPACE  SYSTEMS  ACTIVITY 
P.o.  BOX  92960 
wORLOwAY  POSTAL  center 
lOS  ANGELES,  ca.  90009 
01CY  ATTN  CODE  52 

OFFICE  OF  NAVAl  RESEARCH 
ARLINGTON,  VA  22217 

01CY  ATTN  CODE  465 
01CY  ATT\  CODE  461 
01CY  ATTN  CODE  402 
01CY  ATTN  COOE  420 
01CY  ATTN  CODE  421 

COMMAfOER 

AEROSPACE  DEFENSE  COMMAfD/X 
DEPARTMENT  of  The  air  FORCE 
ENT  AFB,  CO  80912 

01CY  ATTN  DC  MR.  LONG 


COMMANDER 

U.S.  ARMY  MATERIEL  DEV  &  READINESS  CMD 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DRCLX  J.  A.  BENDER 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND /XPD 
DEPARTMENT  OF  T»-E  AIR  FORCE 
ENT  AFB,  CO  80912 

01CY  ATTN  XPDQQ 
QIC Y  ATTN  XP 


COMMANDER 

U.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  BACKLICX  ROAD 
BLDG  2073 

SPRINGFIELD,  VA  22150 
01CY  ATTN  LIBRARY 

DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH  LABS 
ABERDEEN  PROVING  GROUND,  MD  21005 

01CY  ATTN  TECH  LIB  EDWARD  BA ICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MOMOUTH,  NU  07703 

01CY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  Al  35809 
01CY  ATTN  JIM  GAMBLE 


DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS  ACTIVITY 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
01CY  ATTN  ATAA-SA 
01CY  ATTN  TCC/F.  PAYAN  JR. 

01CY  ATTN  ATAA-TAC  LTC  J.  HESSE 


AIR  FORCE  GEOPHYSICS  LABORATORY 

HANSCOM  AFB,  MA  01731 

01CY  ATTN  OPR  HAROLD  GARDNER 
01CY  ATTN  OPR-1  JAMES  C.  ULwIC* 

01CY  ATTN  LK8  *E**€Th  S.  W.  CHAMPION 
01CY  ATTN  OPR  ALVA  T.  STAIR 
01CY  ATTN  PhP  JULES  AARONS 
01CY  ATTN  PK>  JURGEN  BUCHAU 
01CY  ATTN  Pt-D  JOHV  P.  MULLEN 

AF  WEAPONS  LABORATORY 

XIRTLAND  AFB,  nh  87117 
01CY  ATTN  SUL 

01CY  ATTN  CA  ARTHUR  H.  GUENTHER 
01CY  ATTN  DYC  CAPT  J.  BARRY 
01CY  ATTN  DYC  JOHN  ,M.  <AMM 
01CY  ATTN  DYT  CAPT  .MARK  A.  FRY 
01CY  ATTN  DES  MAJ  GARY  GANONG 
01CY  ATTN  DYC  J.  JANNI 

AFTAC 

PATRICK  AFB,  Fl  32925 

01CY  ATTN  TF/MAJ  WIlEY 
01CY  ATTN  TN 


AIR  FORCE  AVIONICS  LABORATORY 
WR1GHT-PATTERSCN  AFB,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
01CT  ATTN  AAD  ALLEN  JOHNSON 


COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 

WASHINGTON,  D.C.  20360 

01CY  ATTN  NAVALEX  034  T.  HUGHES 
01CY  ATTN  PME  117 
01CY  ATTN  PME  117 — T 
01CY  ATTN  COOE  5011 

COmANDING  OFFICER 

NAVAl  INTELLIGENCE  SUPPORT  CTR 

4301  suitland  road,  bldg.  5 

WASHINGTON,  D.C.  20390 

01CY  ATTN  MR.  DUBBIN  STIC  12 

01CY  ATTN  NlSC-50 

01CY  ATTN  COOE  5404  J.  GALET 


DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  (  ACQ 
DEPARTMENT  OF  The  air  FORCE 
WASHINGTON,  D.C.  20330 
01CY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/XR 
DEPARTMENT  of  The  air  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  XR  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DI Vi S I  ON/ Y SEA 
DEPARTMENT  OF  the  A|R  force 
HANSCOM  AF0,  MA  01 731 
01CY  ATTN  YSEA 
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Gc*€RAL  ELECTRIC  TECH  SERVICES  CO.,  InC. 
hMES 

COURT  STREET 
SYRACUSE,  NT  13201 

01CY  ATTN  G.  HILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

01CY  ATTN  JO**  ISE  JR 
01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBAMLS,  AX  99701 

CALL  CLASS  ATTn:  SECURITY  OFFICER) 
01CY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

01CY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  0219** 

01CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
JRBANA,  IL  61803 

01CY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFHARD 
OICY  ATTN  JOEL  SENGS TON 

HSS,  INC. 

2  Alfred  circle 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

JNTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMINO  DEL  MAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JO**S  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 

oicy  attn  Thomas  potew*a 

OICY  ATTN  JOHN  DASSOULAS 

LOCK  FEED  MISSILES  l  SPACE  CO  INC 
P.  0.  BOX  504 
SU**YVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AN©  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  W.  L.  imOF  DEPT  52-12 


KAMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

LINCABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

LOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

OICY  ATTN  K.  BIBL 

M.I.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M,  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 

MARTIN  .MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  80X  5837 
ORLAhOO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

MCDOWELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.  W.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
OICY  ATTN  W.  F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  0.  SAPPENFIELD 
OICY  ATTN  R.  BOGUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCHEI8E 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JO**  MORGANS TERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RF I  ELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.  C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADORESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 


rHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

physical  dynamics  inc. 

P.  O.  0OX  3027 
3ELLEVUE,  wA  98009 

01CY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DY'.AMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

01 CY  ATTN  A.  THOMPSON 

RID  ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  3.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  lELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  .MAIN  STREET 
SANTA  MONICA,  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  30X  2351 
LA  JOLLA,  CA  92038 

01CT  ATTN  LEWIS  M.  LINSGN 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  0.  SACHS 
OiCY  ATTN  E.  A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  ma  01776 

OICY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTSVILLE  DIVISION 
2109  w.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  Al  35805 

OICY  ATTN  DALE  H.  DIVIS 


SCIENCE  APPLICATIONS,  INCORPORATED 
8400  wESTPARK  DRIVE 
MCwEAN,  VA  22101 

OICY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
Pleasanton,  ca  94566 
OICY  ATTN  SZ 


SRI  INTERNATIONAL 
333  RAVENS *000  AvENUE 
MEM.C  PARK,  CA  94025 

OICY  ATTN  DONALD  NEIL SON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.  L.  COBB 
OICY  ATTN  DAVID  A.  oQHNS ON 
OICY  ATTN  WALTER  G.  ChESNUT 

oicy  attn  Charles  l.  ring 

OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  3AR0N 
OICY  ATTN  RAY  L.  LEAD AS RAND 
OICY  ATTN  G.  CAPP-ENTER 
OICY  ATTN  G.  PRiCE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  international  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  017  3C 

OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  t  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONX  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLEBUCH 
oicy  attn  s.  Altschuler 
oicy  attn  d.  dee 

VISIDYNE,  INC. 

19  third  avenue 

NORTfi  WEST  INDUSTRIAL  PARK 
Burlington,  ma  01803 

OICY  ATTN  CHARLE:  HUMPHREY 

oicy  attn  j.  w.  carpenter 
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